scalar maxDi = pyrolysis->maxDiff();
